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abstract  .vC.  .; 

The  acceleration  of  the  mean  lunar  longitude  has  a  small  effect  on  the  periods  of 
most  terms  in  a  Fourier  expansion  of  the  longitude.  There  are  several  planetary 
perturbation  terms  that  have  small  amplitudes,  but  whose  periods  are  close  to  the 
resonant  period  of  the  lunar  libration  in  longitude.  Some  of  these  terms  are  moving 
toward  resonance,  some  are  moving  away  from  resonance,  and  the  periods  of  those  terms 
that  do  not  include  the  Delaunay  variables  in  their  arguments  are  not  moving. 

Because  of  its  acceleration  of  longitude,  the  Moon  is  receding  from  the  Earth,  so  the 
magnitude  of  the  restoring  torque  that  the  Earth  exerts  on  the  rotating  Moon  is 
gradually  attenuating;  thus  resonance  itself  is  moving,  but  at  a  much  slower  rate 
than  the  periods  of  the  accelerating  planetary  perturbations.  There  are  five 
planetary  perturbation  terms  from  the  ELP-200()  Ephemeris  (with  amplitudes  of  0 
.00001  or  greater)  that  have  passed  through  resonance  in  the  past  two  million  years. 
One  of  them  is  of  special  interest  because  it  appears  to  be  the  excitation  source 
of  a  supposed  free  libration  in  longitude  that  has  been  detected  by  the  lunar  laser 
ranging  experiment.  The  amplitude  of  the  term  is  only  0"  .00021  but  it  could  be  the 
source  of  the  approximately  1"  amplitude  free  libration  term  if  the  viscoelastic 
properties  of  the  Moon  are  similar  to  those  of  the  Earth. 


14.  SUBJECT  TERMS 

Lunar  librations,  lunar  dissipation,  lunar  Q,  free  llbrations 

1 

15.  NUMBER  OF  PAGES 

18  ■ 

16.  PRICE  CODE 

,  17.  security  classification 

18.  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATION 

20.  LIMITATION  OF  ABSTRACT  , 

OF  REPORT 

Of  THIS  PAGc 

Of  ABSTRACT 

1  UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

SAR  ; 

Best 

Available 

Copy 


PL-TR-93-2223 


PASSING  THROUGH  RESONANCE:  THE  EXCITATION  AND 
DISSIPATION  OF  THE  LUNAR  FREE  LIBRATION  IN  LONGITUDE 

DONALD  H.  ECKHARDT 
Geophysics  Directorate,  Phillips  Laboratory 
Hanscom  AFB,  MA  01731,  USA 

Abstract.  The  acceleration  of  the  mean  lunar  longitude  has  a  small  effect  on  the  periods  of  most  terms 
in  a  Fourier  expansion  of  the  longitude.  There  are  several  planetary  perturbation  terms  that  have  small 
amplitudes,  but  whose  periods  are  close  to  the  resonant  period  of  the  lunar  libration  in  longitude.  Some 
of  these  terms  are  moving  toward  resonance,  some  are  moving  away  from  resonance,  and  the  periods 
of  those  terms  that  do  not  include  the  Delaunay  variables  in  their  arguments  are  not  moving.  Because 
of  its  acceleration  of  longitude,  the  Moon  is  receding  from  the  Earth,  so  the  magnitude  of  the  restoring 
torque  that  the  Earth  exerts  on  the  rotating  Moon  is  gradually  attenuating;  thus  resonance  itself  is 
moving,  but  at  a  much  slower  rate  than  the  periods  of  the  accelerating  planetary  perturbations.  There 
are  five  planetary  perturbation  terms  from  the  ELP-2000  Ephemeris  (with  amplitudes  of  0"  .00001  or 
greater)  that  have  passed  through  resonance  in  the  past  two  million  years.  One  of  them  is  of  special 
interest  because  it  appears  to  be  the  excitation  source  of  a  supposed  free  libration  in  longitude  that 
has  been  detected  by  the  lunar  laser  ranging  experiment.  The  amplitude  of  the  term  is  only  0"  .0(X)21 
but  it  could  be  the  source  of  the  ~  1 "  amplitude  free  libration  term  if  the  viscoelastic  properties  of 
the  Moon  are  similar  to  those  of  the  Earth. 

Key  words:  Lunar  librations,  lunar  dissipation,  lunar  Q.  free  librations 


1.  The  Rigidity  of  the  Moon 

The  velocity  of  a  seismic  shear  wave  is 


where  p  is  the  rigidity  and  a  is  the  density  of  the  medium.  The  structure  of  the 
lunar  interior  is  rather  homogeneous.  Over  the  depth  range  60-1 100  km,  the  shear 
wave  velocity  varies  little  from  Vs  =4.3  km/s  in  the  model  (designated  GDT)  of 
Goins,  Dainty  and  Toksoz  [1981];  Vg  is  about  five  per  cent  higher  over  the  same 
depth  range  in  the  model  (designated  N)  of  Nakamura  [1983].  Let  M  be  the  mass 
of  the  Moon,  R  its  radius,  C\\  its  minimum  moment  of  inertia,  C33  its  maximum 
moment  of  inertia,  and  C22  the  moment  of  inertia  about  the  third  principal  axis. 
Because  /?  =  [C33  -  Cn]/C22  =  0.000631  and  C^^/MR^  =  0.392  [Williams, 
1977],  all  the  moments  of  inertia  are  close  to  that  of  a  homogeneous  sphere,  so 
the  lunar  density  {a  =  3340  kg  m“^)  is  fairly  uniform.  Thus,  the  lunar  rigidity, 
calculated  using  (1),  does  not  vary  much  from  p  =  62  x  lO’  Pa  (GDT  model)  or 
/li  =  68  X  10^  Pa  (N  model).  (The  lunar  rigidity  is  somewhere  between  that  of  cast 
iron,  p  =  40  X  10^  Pa,  and  steel,  p  =  iOx  10’  Pa.) 

For  an  incompressible  non-dissipati ve  homogeneous  spherical  Moon,  its  second 
degree  potential-disturbance  Love  number  (omitting  the  subscript  in  ^2)  is  k  = 
3/[2  +  \9p/agR]  [Love,  pp.  249-259, 1927;  Jeffreys,  p.  299, 1976],  where  g  is  the 
acceleration  of  gravity  at  the  lunar  surface.  Using  agR  =  9.5  x  lO’  Pa,  the  Love 
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number  determined  from  seismic  data  is  k  =  0.024  (GDT  model)  or  k  =  0.022 
(N  model).  Cheng  and  Toksoz  [1978]  used  numerical  integrations  to  calculate  the 
Love  number  for  two  compressible  Moon  models  in  which  the  elastic  properties 
are  functions  of  depth.  The  differences  between  their  estimates,  k  =  0.029  and 
k  =  0.034,  and  those  from  the  GDT  and  N  models  are  principally  due  their  use  of 
earlier  shear  velocities  which  are  appreciably  lower  than  the  more  accurate  GDT 
and  N  velocities. 

The  Moon  deforms  elastically,  causing  the  inclination  of  its  equator  to  the 
ecliptic  to  increase  by  61  =  kx  74"  [Eckh^t,  1981].  From  analyses  of  lunar  laser 
ranging  data,  Williams,  Newhall  and  Dfckey  [1987]  estimate  that  k  =  0.027  ± 
0.006;  that  is,  the  deformation  causes  ^e  inclination  to  increase  by  2"  resulting 
in  a  -17  m  sinF  term  (7’  is  the  argument  of  latitude)  in  the  latitude  libration. 
Other  libration  parameters  (e.g.,  /)  and  Cn/M R^)  contribute  to  the  sin  F  term 
in  the  latitude  libration,  but  they  contribute  to  other  terms  in  different  ways,  so  a 
thorough  analysis  allows  the  deformation  term  to  be  discriminated. 

To  take  into  account  anelastic  dissipation  in  the  Moon,  the  Love  number  is 
replaced  in  the  frequency  domain  by  a  complex  Love  number,  k{  1  -  i/Q ),  where 
Q  >  1 .  The  term  {  \  -  i/Q)  a  exp( —i/Q)is  actually  a  phase  shi  fting  operator  that 
changes  a  periodic  term  of  the  form  exp(  iujt)  to  exp[i(a;t  -  1  /Q )].  The  imaginary 
part  of  the  Love  numbci  is  negative  because  a  deformation  lags  its  indueing  forces. 
The  magnitude  of  the  ratio  of  the  real  part  to  the  imaginary  part  of  ^(1  -  i/Q) 
is  Q,  where  l/<5  is  the  specific  dissipation  function  [Munk  and  MacDonald,  pp. 
21-22, 1960]  which  is  frequency  dependent.  Because  \9fi/agR:$>  2,  the  following 
approximation  is  valid  for  the  Moon 


19/iA;  S3  3agR. 


(2) 


CD 

CO 


Because  >  1 ,  (2)  remains  valid  if  A’  is  replaced  by  k{]  -  i/Q)  and  /r  is  replaced 
by  /i(l  -|-  i/Q).  The  imaginary  part  of  the  rigidity  is  positive  because  the  stress 
leads  the  strain  [Ben-Menahem  and  Singh,  pp.  856-859, 1981  ].  A  value  for  k  will 
subsequently  be  required  in  order  to  estimate  Q]  taking  into  account  the  various 
estimates  discussed  above,  this  paper  adopts  the  value  k  =  0.025. 


2.  Equations  Describing  the  Libration  in  Longitude  of  a  Deformable  Moon 

In  a  selenocentric  coordinate  system,  let  r;  =  rw,  locale  the  Earth’s  center  of  mass, 
a  distance  r  from  the  origin,  and  let  pi  =  pvi  locate  a  point  within  the  Moon,  a 
distance  p  from  the  origin.  The  lunar  gravitational  potential  at  r,  is 

=  {Glr)[  f  dM  +  f^iR/r)^  f  {p/ R)^Pn(u,v,)dM] 

JMoon  n=2  JMoon 


where  G  is  the  gravitational  constant,  P„  is  the  Legendre  function  of  degree  n,  and 
the  integrals  are  over  the  mass  of  the  Moon.  Through  degree  n  =  2,  the  potential 

93-27948 


LUNAR  FREL  LIBRATLON  IN  LONGITUDE 


309 


is 

V  =  GMtr  +  {Gl?)[\  I  -  p,'^]dM  +  /  p,^dM]. 

^  JMoon  JMoon 

Using  the  moment  of  inertia  tensor, 

C{j  —  /  (PA:  piPi)dM ^ 

J  M  oon 

where  6ij  is  the  Kronecker  delta,  the  truncated  potential  is  reformulated  as 

V  =  GM/r  -  {G /2r^){3u,C,jUj  -  C,.)- 

The  potential  and  moments  of  inertia  of  an  anelastic  Moon  are  deformed  by  forces 
derived  from  its  centrifugal  potential,  Wi,  and  from  the  second-degree  Earth- 
induced  tidal  potential,  W2.  The  lunar  potential  at  the  Earth  is  perturbed  by 

6V  =  lc(\-  ilQ)(Rlr)\Wx  +  W2} 

where 

1,1  , 

ITi  =  -  ^ij<^k  )rj, 

and 

W2  =  {Gm/2r^)r,{3uiUj  -  Sij)rj, 

where  Ui  is  the  lunar  angular  rotation  rate,  m.  is  the  mass  of  the  Earth,  and  Sijk  is 
the  alternating  tensor.  Then 

6V  =  -{GI2r^){3u,SC\jUj  -  6Ckk^,j)  -1-  k{\  -  i/Q)RW/3r^ 

where  the  perturbation  of  the  moment  of  inertia  tensor  is 

6Cij  =  (\ -  i/Q)K[uiUj  -  3(Gm/a^)(a/rfiiiUj] 

=  (\ ~  i/Q)K[u;iUjj  -  3n^{a/r)^UiUj]  (3) 

%  (1  -  i/Q)K{uj,Ljj  -  3n^UiUj), 


a  is  the  mean  Earth-Moon  distance. 


K  =  kRy3G  = 


k(MR^)iRlaf 

{M/m){3ri^) 


(4) 


and  n  is  the  mean  rate  of  lunar  motion. 

For  a  small  lunar  rotation,  Sipk,  the  Earth’s  direction  cosines  change  by  = 
£ijkUjSipk,  so  dui/dfpk  —  SijkUj.  The,  torque  exerted  on  the  Moon  is 


Nk  =  md{V  -I-  6V)fdipk  =  (3Gmlr^)ui{Cij  +  fiCij)£jiktii.  (5) 
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In  an  inertial  frame,  the  time  derivative  of  the  angular  momentum,  L,  -|-  6L,  = 
{Cij  +  6Cij)ujj,  is  equal  to  the  torque,  so 

d{Li  +  iLi)j dt  =  Li  -]r  f>Li  +  Sijk*^j(Lk  +  ^Lk)  =  Ni,  (6) 

where  the  dot  derivative  is  with  respect  to  time  in  the  lunar  coordinate  system 
which  is  chosen  to  be  the  same  as  the  principal  axes  of  Cij.  Setting  (5)  equal  to  (6) 
with  SCij  =  0  gives  the  basic  Euler  dynamical  equations  for  the  rotation  of  a  rigid 
Moon. 

If  the  Love  number  is  real,  the  centrifugal  potential  deforms  the  Moon  as  an 
oblate  spheroid  whose  shortest  axis  is  the  instantaneous  rotation  axis,  so  \V\  docs 
not  affect  the  rotation  rate  through  the  RHS  of  (6),  only  its  direction.  Also  if  the 
Love  number  is  real,  the  tidal  potential  deforms  the  Moon  as  a  prolate  spheroid 
with  its  longest  axis  in  the  instantaneous  direction  of  the  Earth,  so  there  can  be 
no  torque  due  to  W2.  The  real  part  of  the  Love  number  only  has  an  effect  on  the 
lunar  rotation  rate  through  the  change  in  the  Moon’s  moment  of  inertia,  and  that  is 
almost  entirely  due  to  the  centrifugal  potential.  The  mean  lunar  rotation  rate  n  (the 
same  as  its  mean  rate  of  motion)  is  much  larger  than  its  perturbation  r,  the  libration 
in  longitude,  so  the  real  part  of  the  Love  number  affects  the  Moon’s  librations 
principally  through  a  “frozen  in’’  deformation  in  its  moment  of  inertia  tensor  which 
is  a  permanent  part  of  Cij.  The  impact  of  the  real  part  of  the  Love  number  on 
the  rotation  rate  is  essentially  nil;  only  the  imaginary  part  contributes  significantly 
to  perturbations  from  the  Euler  dynamical  equations.  This  can  be  demonstrated 
mathematically  by  combining  (3),  (4),  (5)  and  (6). 

For  examining  the  libration  in  longitude  r  (about  the  3-axis),  use  the  approxi¬ 
mations 


■  0  ■ 

■  1  ■ 

OJi  = 

0 

,Ui  = 

U2 

.‘^3  . 

.  0  . 

UJ  =  n  f  &  U2  —  S  —  T, 

where  the  1-axis  is  in  the  mean  Earth  direction  (or  nearly  so  [Eckhardt,  19731)  and 
s  is  comprised  of  the  equation  of  the  center  and  all  inequalities  in  longitude.  With 
7  =  (<^22  -  C'h)/C'33  =  0.(X)0228  and,  approximately,  the  bottom 

row  of  the  Euler  equation  becomes  approximately 

f  =  u>o^(s  -  t).  (7) 

Excitation  and  dissipation  require  hysteresis,  so  orJy  the  imaginary  periodic  terms 
of  the  moment  of  inertia  perturbation, 

— 3n^  — 3n^W2  0 

—3n^U2  —3n^U2^  0  , 

0  0 


sc,j  =  {I  ~i/Q)ri 
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are  significant.  Therefore,  to  the  first  order  in  s,  r  and  f,  the  significant  portion  of 


SCij  is 


=  -iiK/Q) 


0  -3n^(s  —  t)  0 

—3n^(s  -  r)  0  0 

0  0  2nf 


Dropping  second  order  terms,  the  variable  perturbation  in  the  lunar  angular  mo¬ 
mentum  is 


=  SCijUj  —  -2i{KIQ)n?’f 
and  its  rate  of  change  is 

SLi  =  SCijUj  =  —2i{K  IQ)Tt^f 

Then 

iCijUj  =  3i{KIQ)n^ 

and 

SkiitLiSC’ijUj  =  3i(I\/Q)n^ 


0 

S  —  T 

0 


0 

0 

s  —  T 


The  perturbation  of  the  torque  exerted  on  the  deformed  Moon  by  the  Earth  is  then 

0 

6Ni  =  m€ijkUj[d{6V)/duk]  =  9i{K/Q)n*  0 

S  -  T 


The  perturbation  of  (7)  accounting  for  the  deformable  Moon  model  is 

{I -2i'yK/3)f  =  {\  +  iK)oJo^{s  -  t)  (8) 

where  k  =  <;k/Q  and  t  =  (Rfa)^{MR?-fC$3)‘y~^{m/M)  =  0.0841.  Take  k  = 
0.025  and  suppose  that  Q  >  20  so  that  k  <  lO"**.  Because  27  C  3,  the  2i7k/3 
term  is  dropped,  and  (8)  becomes 

f  =  a;o^(l -I- ?;k)(5  -  r).  (9) 

The  period  of  resonance,  27r/a;o,  derived  by  detailed  semi-analytic  theories  has  a 
small  uncertainty.  Estimates  based  on  the  LURE  2  constants  [Williams,  1977]  are. 
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in  sidereal  months,  38.651  [Eckhardt,  1981],  38.654  [Moons,  1982]  ,  and  38.658 
[Migus,  1980].  These  estimates  are  adjusted  for  the  JPL  constants  of  Williams, 
Newhall  and  Dickey  [1987].  The  increase  in  7  is  (227.951  -  227.370)  x  10“*  = 
0.581  X  10“*  which  decreases  each  estimate  by  0.049,  but  using  the  technique  of 
Williams  et  al.  [  1973],  the  changes  in  the  harmonic  coefficients  C31  and  C33  work 
out  to  be  the  same  as  decreasing  7  by  0.592  x  10“^  and  the  adjustments  almost 
cancel  each  other.  The  net  result  is  that  each  estimate  increases  by  only  0.001 ;  this 
reflects  the  fact  that  the  lunar  laser  ranging  experiment  is  more  sensitive  to  the 
period  of  r  resonance  than  to  the  related  solution  parameters. 

3.  Analytic  Solutions  of  the  Libration  Equations 
Let  the  longitude  source  terms  be 

where  each  //*,.  is  a  complex  coefficient.  Then  the  solutio  (9)  is 

_  ( 1  +  iK)Sk  _ Sk _ 

^  ^  \  +  iK-  (w/t/wo)^  ^  1  -  (u)k/u)o)'^  ‘ 

The  i-K  term  has  negligible  effect  on  the  steady  state  solution  given  by  (10)  because 
currently  none  of  the  Sk  source  terms  has  a  frequency  close  enough  to  resonance. 
If,  however,  the  argument  of  some  source  term  includes  an  acceleration  term, 

Sk  =  Hke\p[i(ukt  +  (1 1) 

the  dissipation  term  r.iay  be  important,  especially  near  resonance.  Perturbing  this 
process  is  the  fact  that  because  uq  is  proportional  to  n,  it  is  a  linear  function  of 
time  (because  of  the  acceleration  of  the  mean  longitude  of  the  Moon). 

Suppose  that  there  is  a  term  of  the  form  (11)  that  carried  t  through  resonance 
in  the  recent  past  (using  a  geological  sense  of  time  so  that  a  million  years  or  so 
ago  was  recent)  at  <  =  0  so  that  u^iO)  =  wo(0).  Change  the  time  variable  to  the 
dimensionless  variable  77  where  fj  =  0^0(0)  so  that  rj  goes  through  27r  radians  every 
38.65  sidereal  months,  and  use  the  prime  (')  to  denote  the  derivative  with  respect 
to  77.  The  argument  of  s/t  in  (1 1)  is  then  expressed  in  the  form 

+  ^Ckt^  =  u;o(0)t  +  ^Cki^  =  77  + 
in  which  case 

Xk  =  CklijJQ^iO). 

In  the  Fourier  expansion  of  the  longitude  of  the  Moon,  many  of  the  planetary  source 
terms  near  resonance  have  accelerations  with  approximately  the  same  magnitude 
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as  the  acceleration  of  the  mean  longitude  of  the  Moon  [Williams,  Newhall  and 
Dickey,  1992], 

n  =  (-26.0  ±  1.0)7cy 

or 

h/n  =  -0.150  X  10"’/y. 

These  are  the  only  plarxtaiy  terms  of  any  signihcance  for  this  study.  Thus,  with 
u;o(0)  =  2.17/y, 

\Xk\  «  26".0cy"V‘^^(0)  =  2.67  x  10"^  (12) 

Let  the  resonance  angular  rate  as  a  function  of  r/  be  represented  by  with  this 
choice,  fc'(O)  =  1.  Also  define  Xq  =  l>'.  Because  i>  is  proportional  to  n, 

h  dlnn(f)  _  dlnvin) 
and,  therefore, 

Ao  =  -0.150  X  10-7"V^o(0)  =  -0.069  x  10"^  (13) 

Then,  on  dropping  the  negligible  term, 

°  =  a;o^(0)[l  +  “  a;o^{0)(l  +  2A077). 

Change  the  independent  variable  of  (9)  from  t  to  r}, 

r'k  +  (1  +2Ao??)(l  +  in)Tk  =  (1  +2Ao7?)(l  +  iK)HkCxp[l{r)  +  Aa:?7^/2)].(14) 
Next  substitute 

=  6(7)exp(tr/)  (15) 

into  (14),  where  ^kiv)  'S  a  function  that  varies  slowly  compared  with  exp(i77)  so 
that  |^fc(77)|  <  |^^77)|.Then,  on  neglecting  and  all  terms  quadratic  in  k,  Xq,  Hk, 

^k  +  (k/2  -  iXoi])^k  =  -i{Hkl2)e\p{iXkV^/2).  (16) 

Using  the  integrating  factor  exp[(K7;  -  iXoif)/!],  the  solution  to  (16)  at  r/  =  /i  is 

C  *1/  j-Kh  +  iXoh'^  ri{Xk  -  Xq)t]^  +  KTf 

^kih)  =  --Hkexpi - - - ]  J  exp[-i - - - Jdr/.  (17) 

From  (12)  and  (13),  it  is  clear  that  A^  dominates  the  A^.  -  Aq  term.  The  integral  for 
Ao  >  Xk  (Xk  obviously  being  negative)  is 

I  exp[-^^^5 — =  exp[-f(K/pfc)^/8]  /  exp(-u;^)d7?  (18) 
J—oo  ^  J—oo 
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where 

w  =  w{r])  =  (1  +  i)Pk‘nl'^  -  (1  -  i)f^l{^Pk)  (19) 

and 

Pit  =  |Ao-A,|'/^  (20) 

(The  absolute  value  is  used  in  this  definition  so  that  it  will  apply  as  well  when  At 
is  positive.)  Substituting  (1  -  i)pk~^dw  =  dt], 

/h.  1—1  ^  1—1 

c\p(-'uP')dri  — -  /  exp(-ti;^)dtw  =  — — \/^[l  +  erf[w(/i)]],  (21) 

-oo  Pit  J —oo  ^Pk 

where  erf[u;(/i)]  is  the  error  function  with  the  complex  argument  iv(h). 

The  phase  of  Hk  is  the  phase  of  Sk  precisely  at  t  =  0.  There  is  no  way  of 
knowing  the  phase  of  Hk,  so  there  is  no  way  of  knowing  the  phase  of  Tk  near 
resonance.  Only  the  magnitude  of  Tk  can  be  calculated  as  a  function  of  \  Hk\,  Pk,  « 
and  h\  combining  (15),  (17),  (18)  and  (21),  it  is 

|n.(/i)|  =  yV2[|/ffc|/(2p,)]exp(-K/i/2)|l  +crf[uK/i)]|.  (22) 

The  integral  (1 8)  for  <  Aq  is  the  complex  conjugate  of  the  integral  for  A*..  >  Aq, 
so  (22)  is  valid  for  all  values  of  Aq  -  A*,  except,  of  course,  for  (the  non-existent) 
\k  =  Ao  which  represents  being  parked  at  resonance. 

Hlpk^h  >  K, 

w{h)  ^  i)Pkhl2  (23) 

and 


^^V^[l  +erf[w(/i)]]  w  ^[[\  +2C(p/t/»/\A)]  -  »[1  +  2S(pfc/i/>A)]], 

2pk  2pk 

where  C{pkh/ y/n)  and  S{pkh/ x/x)  are  the  Fresnel  integrals.  Then 

|rt(/i)|  w  V^['^]®xp(-y)y/[l  +2C(^)P  4-  [1  +  2S(^)]2.  (24) 

For  pkh  >  y/ir,  C^pkh/y/n)  «  S{pkh/ y/ir)  w  1  /2,  so  (24)  becomes 

|rfc(/i)|  «  yfT^[\Hk\/pk]c\p(-Kh/2).  (25) 
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4.  The  Accelerating  Arguments 

Let  m  —  n'/n  =  0.0748  currently  where  n'  is  the  mean  rate  of  solar  motion. 
Numerical  integrations  of  the  orbits  of  the  Solar  System  planets  demonstrate  that 
their  mean  motions  have  inappreciable  accelerations  over  millions  of  years  (Quinn, 
Tremaine  and  Duncan,  1991],  so  n'  =  0  and  rh/m  =  —njn.  Using  the  m  power 
series  for  the  mean  motions  of  the  lunar  node  Q  and  perigee  u  (Brouwer  and 
Qemence,  pp.  322-323, 1961],  their  accelerations  are  then  related  to  h  by 

Q  =  (1.5m2  -  0.84375m^  -  8.531252m''  -  23.91 845m5 
-48.6504m^  -  79.016m'')n  =  0.0077n 

and 

=  (-1.5;n2  -  21.0312571)3  -  127.218752m''  -  648.17625m3 
-3130.5256m^  -  15118.8877/)''  -  80282.32/))* 

-42578.37))^)/)  =  -0.0236/). 

The  accelerations  of  all  components  of  the  planetary  perturbation  arguments  arc 
essentially  zero  except  for  the  Delaunay  variables  D,  I  and  F  which  have  the 
accelerations,  respectively,  h,h  -  Cj  =  1.0236/),  and  h  -  Hi  =:  0.9923/);  in  units 
of  //,  they  are  h/u!Q^{0),  1 .0236r)/u)o^(0),  and  0.9923/)/u)o^(0)  where  /)/u;o^(0)  = 
-2.67  X  10"^.  For  a  Fourier  term  in  the  planetary  longitude  perturbation  of  the 
Moon  to  be  of  interest,  its  period  must  be  near  resonance  and,  with  its  argument 
expressed  as  an  increasing  function  of  time,  the  sum  of  the  integer  coefficients 
of  D,  I  and  F  in  the  argument  must  be  equal  to  (or  greater  than)  1  if  the  term’s 
period  is  greater  than  the  resonant  period  or  the  sum  must  be  equal  to  (or  less 
than)  -1  if  its  period  is  less  than  the  resonant  period.  The  effect  of  the  slowly 
changing  resonance  (the  Aq  term)  can  be  overlooked  in  this  winnowing  compared 
with  the  uncertainty  in  the  period  of  resonance,  a  value  that  is  theory  and  parameter 
dependent.  Source  terms  from  the  Chapront  and  Chapront-Touzc  (1983]  ELP-20(X) 
lunar  ephemeris  that  meet  these  criteria  are  given  in  Table  \.V  ,T  and  M  represent 
the  mean  longitudes  of  Venus,  Earth  and  Mars. 

The  nominal  period  of  resonance  is  chosen  as  38.652  sidereal  months  (Eckhardt, 
1981];  adjustments  for  the  resonant  periods  of  Migus  (1980]  and  Moons  (1982] 
will  be  considered  later.  Because  of  their  amplitudes  and  closeness  to  resonance 
(one  above  and  the  other  below),  the  last  two  terms  in  Table  1  are  of  principal 
interest .  The  other  three  are  insignificant  and  will  be  neglected.  For  the  subscripts 
designating  these  terms,  use  1  for  16331  and  2  for  16018. 

Term  16331  is  closest  to  resonance  and,  by  far,  the  most  important  of  all 
terms  that  have  “recently”  passed  through  resonance.  Its  acceleration  is  (2  - 
1.0236)f)/u;o^(0),  so 

Ai  =  0.9764  X  (-2.67  x  10“^)  =  -2.61  x  10"’, 
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TABLE  I 

Planetary  perturbations  in  longitude  for  terms  that  have  “recently”  passed  through  resonance 
and  whose  periods  differ  from  the  resonance  period  by  less  than  one  per  cent.  These  terms 
were  selected  from  the  ELP-2000  lunar  ephemeris  of  Chapront  and  Chapront-Touz6  [1983]. 
The  sequence  numers  are  those  of  ELP-2000 


Sequence 

Number 

m 

(Arc-Seconds) 

Argument 
(Excluding  Phase) 

Period 

(Sidereal  Months) 

0.00011 

-39V  +  39r  -f  2D 

39.080 

16766 

0.00003 

-24  V  +  30T  -4M  +2D-I 

39.039 

16979 

0.00001 

-26V  +  29T  +  4M  +2D-1 

16331 

0.00021 

-21V  +23T-^2D-l 

38.666 

16018 

0.00066 

19  V  -  18T  -  D -  F 

38.258 

and  the  Aq  term  just  marginally  contributes  to  making  [see  (13)  and  (20)] 

=  2.54  X  10'^ 

Then,  because  the  dimensionless  angular  rate  at  resonance  is  unity, 

1  -f  pi^/i  =  38.666/38.652  =  1  -t-  0.000362 
and 

h  =  143000. 

Resonance  occurred  (2.8912/27r)/i  =  66(X)0  years  ago.  At  that  time,  its  libration 
amplitude  was  half  of 

y^^lff,|/pi=5".2 

and,  assuming  that  (25)  affords  a  close  enough  approximation,  its  amplitude  is  now 

|r,(/i)|  w  5".2exp(-71000K).  (26) 

The  more  precise  equation  (22)  will  be  applied  later.  Following  the  same  calcula¬ 
tions  for  Term  16018  gives 

\T2ih)\  K  16'M  exp(-1.9  X  IO^k).  (27) 

From  the  ratio  (26)/(27), 

In(/i)l/k2(/i)|  >  1  for  K  >  0.62  X  10"^.  (28) 

From  the  analysis  of  a  six  year  span  of  lunar  laser  ranging  data,  Calame  [1977] 
estimated  the  free  libration  in  longitude  to  have  a  1".8  amplitude.  She  did  not. 
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Fig.  I .  Plots  of  I  n  1  at  the  current  epoch  as  functions  of  k.  For  the  Eckhardt,  Moons  and  M  igus 
resonances,  the  solutions  to  In  |  =  1  are  indicated  at  k  =  23.2  x  10~*,  k  =  34.6  x  10“^, 
and  N  =  40.1  x  10“*. 

however,  allow  for  the  planetary  terms  in  lunar  longitude.  Except  for  T\(h)  and 
T2{h),  the  neglected  terms  (principally  ELP-2(XX) Term  6649  which  has  amplitude 
0".85,  argument  3T  -  5M  and  period  39.138  sidereal  months  [Eckhardt,  1982]) 
account  for  only  about  half  of  her  1".8  estimate.  Suppose  that  the  approximately 
1 "  .0  residual  is  due  to  rj  ( /i )  or  r2(  /i ).  If  it  were  entirely  due  to  T\(h)  then,  by  (26), 

K  =  23  X  10-^  (29) 

and  if  it  were  entirely  due  to  T2(h)  then,  by  (27),  k  =  1.5  x  10“^;  but,  by  (28), 
|ri(/i)|  >  lr2(/))j  if  K  =  1.5  X  10~^,  so  (29)  provides  the  single  estimate  of  the 
dissipation  parameter  for  the  Eckhardt  [1981]  resonance  and  1".0  residual.  This  is 
only  an  approximate  solution  and  modifications  are  also  required  for  the  Moons 
[1982]  and  Migus  [1980]  resonances.  Moreover  the  1".0  residual  is  only  a  nominal 
value.  A  plot  of  |ti(143000)|  [Eckhardt],  |ti(112000)|  [Moons],  and  |ri(71000)l 
[Migus]  as  functions  of  k  by  the  general  solution  (22)  is  presented  in  Figure  1 .  The 
solutions  to  |7-i(/*)l  =  1  are  indicated  at 

K  =  23.3  X  10"^  [Eckhardt], 

K  =  34.6  X  10"*  [Moons], 
and 

K  =  40.1  X  10"*  [Migus]. 


(30) 
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Fig.  2.  The  solution  to  (25)  for  ELP-2000  Term  16331  from  132000  years  ago  until 
the  present.  Resonance  according  to  the  adopted  model  occurred  66000  years  ago.  The 
coordinate  units  are  arc-seconds.  There  exists  no  free  libration  with  an  amplitude  as  large 
as  S"  as  shown  in  this  solution,  so  the  effect  of  dissipation  must  be  included  to  get  a 
plausible  model. 

Ah'  Moons’  resonance  is  closer  to  Eckhardt’s  than  Migus’,  her  solution  is 
closer  to  Migus’  than  Eckhardt’s  at  the  current  epoch.  To  understand  how  this  can 
be,  and  to  gain  insight  into  the  nature  of  the  solutions,  the  solution  paths  in  the 
complex  plane  are  examined. 


5.  The  Solution  Paths 

If  K  =  0  and  Ao  =  0,  Equation  (16)  becomes 

^'k  =  -j(//fc/2)exp(tAfc77^/2).  (31) 

Because  |^J.|  =  |//*,/2|,  ^k{i]}  follows  a  path  in  the  complex  plane  at  a  constant 
rate  with  respect  to  -q.  The  direction  of  the  path  is  given  by  the  the  sum  of  the  phase 
of -iHk  and  the  argument  of  the  exponential  term.  Suppose  that  A^  <  0;  then  for 
q  <  0,  the  path  is  clockwise,  and  for  77  >  0,  the  path  is  counter-clockwise.  The 
angular  rate  is  proportional  to  q^  so,  for  r;  <  0,  the  rotation  rate  decreases,  and 
for  q  >  0,  the  rotation  rate  increases.  The  net  result  if  -iHk  is  real  and  positive 
is  that  the  solution  follows  the  Cornu  spiral  [Bom  and  Wolf,  pp.  430-433,  1970] 
depicted  in  Figure  2.  (The  rotations  are  in  the  opposite  directions  for  At  >  0  and 
the  orientation  of  the  spiral  depends  on  the  phase  of  —?//*.) 

If  K  =  0  and  Hk  =  0,  Equation  (16)  becomes 

^k  =  (32) 

In  this  case,  =  |Cfc|  exp(zAo77^/2)  follows  a  circular  path  in  the  complex  plane 
with  an  angular  rate  proportional  to  |  Ao»?i.  The  path  is  counter-clockwise  for  Aq  <  0, 
and  it  is  clockwise  for  Aq  >  0.  A  combination  of  (31)  and  (32)  results  in  a  Cornu 
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Fig.  3.  The  solution  to  (20)  for  ELP- 2000 Term  16331  with  the  same  units  and  over  the  same 
time  span  as  Figure  2.  The  model  is  modified  by  including  a  dissipation  term.  The  current 
1 "  amplitude  of  ri  ( 143000)  shown  in  this  plot  would  be  perceived  as  a  free  libralion. 


spiral  with  the  same  dimensions  as  as  (31);  the  only  effect  of  (32)  is  to  rotate  the 
Cornu  spiral  without  changing  its  dimensions.  For  the  current  problem,  (32)  is  only 
of  marginal  interest. 

If  Ao  =  0  and  Xk  =  0,  Equation  (16)  becomes 

^1.  =  -(k/2)^/,,  (33) 

so  ^kiv)  follows  an  exponential  decay  path  in  the  complex  plane  that  is  a  straight 
line  directed  toward  the  origin.  In  the  general  solution  of  (16),  ^kiv)  t>oth  spirals 
as  in  (31)  and  decays  as  in  (33).  Well  before  passing  through  resonance,  when 
PkV  <  -1  .  ^k{rj)  moves  in  a  tight  spiral  around  the  origin  and  the  decay  term 
has  little  effect.  As  PkP  becomes  small  and  ^kip)  moves  through  the  region  of  the 
resonance  inflection  point,  (31)  dominates  (33).  Eventually,  when  pki]  >  1,  ^kip) 
moves  in  an  increasingly  tight  spiral  around  a  focus  that  decays  toward  the  origin, 
so  (33)  dominates  (31).  The  general  solution  to  (16)  is  given  by  (22)  which  can 
conveniently  be  evaluated  and  plotted  using  Mathematica  [Wolfram,  1988].  A  plot 
of  the  the  T{p)  solution  path  according  to  (22)  (retaining  the  relative  phase  from 
the  error  function)  is  shown  in  Figure  3.  The  dissipation  parameter  used  is  given 
by  (30). 

For  the  shorter  time  intervals  and  higher  dissipations  that  are  demanded  by  the 
Moons  and  Migus  resonances  to  attain  a  1".0  amplitude  for  T\{h),  their  solution 
curves,  shown  in  Figure  4,  are  markedly  different.  The  two  curves  in  Figure  4  .■>re 
currently  headed  in  opposite  directions.  These  directions  arc  relatively  insensiti  ve 
to  either  k  used.  As  time  increases,  the  jri  |  in  the  top  plot  is  decreasing  and  the 
|ti  I  in  the  bottom  plot  is  increasing.  Thus  the  entire  curve  for  the  Moons  solution 
(|ri|  =  l".0atK  =  34.6  X  10“^)  is  moving  downward  with  time  whereas  the  entire 
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Fig.  4.  The  solutions  to  (20)  for  ELP-2000  Term  16331  using  the  Moons  (top)  and  Migus 
(bottom)  resonances  and  k  given  by  (35.2)  and  (35.3).  The  span  of  the  top  plot  is  104(X)0 
years,  and  that  of  the  bottom  plot  is  66000  years. 


curve  for  the  Migus  solution  (|ri  |  =  1".0  at  k  =  40.1  x  10“^)  is  moving  upward 
with  time.  They  only  happen  to  be  close  to  each  other  at  the  present. 

Figure  5  is  a  plot  of  the  libration  amplitude  ]  ri  ( /i)  |  as  a  function  of  the  dissipation 
parameter  k  and  the  resonance  period  in  sidereal  months.  One  can  see  from  this  plot 
that  the  Figure  1  solution  for  the  Eckhardt  resonance  (38.652  sidereal  months)  falls 
in  a  trough;  the  solution  for  the  Moons  resonance  (38.655  sidereal  months)  runs 
along  a  ridge;  and  the  solution  for  the  Migus  resonance  (38.659  sidereal  months) 
falls  in  a  trough.  A  change  of  resonance  of  only  0.001  sidereal  months  moves  a 
solution  from  a  ridge  to  a  trough.  The  resonance  period  is  not  known  well  enough 
to  decide  whether  the  solution  runs  along  a  ridge,  in  a  trough,  or  somewhere  in 
between. 
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Fig.  5.  The  amplitude  of  T\{h)  as  a  function  of  the  resonance  period  (in  sidereal  months) 
and  the  dissipation  parameter,  k. 


Fig.  6.  Plots  of  |ri(/i)|  as  a  function  of  the  resonance  period  (in  sidereal  months)  for  Q 
running  from  30  to  ICiO. 
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The  Steady  state  solution  (10)  for  T2(  h)  has  an  amplitude  of  0".032.  If  approx¬ 
imation  (25)  is  used  to  estimate  |r2(/j)|,  the  result  is  zero,  independent  of  k.  If, 
however,  (22)  is  used,  the  amplitude  is  0".032,  again  independent  of  k.  The  failure 
of  the  approximation  arises  from  the  elimination  of  the  relatively  very  small  term 
( 1  -  i)K/{4pk)  in  replacing  (19)  by  (23).  Because  the  argument  of  the  error  func¬ 
tion  is  not  precisely  7r/4  and  w{h)  is  very  large,  the  error  function  has  a  very  large 
real  part  that  counters  the  very  small  exp(-K/i/2)  attenuation  term.  What  actually 
happens  is  that  ^kiv)  revolves  about  the  origin,  as  indicated  by  the  solution  path  of 
(31),  and  the  revolution  rate  is  just  the  right  amount  to  offset  tJo  to  tJ2.  The  revolu¬ 
tion  velocity  is  slightly  perturbed  by  a  centripetal  velocity  (not  acceleration)  that, 
in  effect,  changes  the  solution  phase  just  as  the  neglected  ik  term  would  change 
the  solution  phase  in  (10).  The  steady  state  solution  is  appropriate  for  ELP-2(KX) 
Term  1 601 8,  so  the  only  term  for  the  libration  in  longitude  that  should  be  modified 
to  allow  for  having  passed  through  resonance  is  the  one  due  to  ELP-2000  Term 
16331.  The  steady  state  libration  solution  for  this  term  should  be  replaced  with  a 
residual  free  libration  term  with  unknown  amplitude  and  phase. 

6.  The  Q  of  the  Moon 

Seismological  and  astronomical  techniques  have  been  used  to  derive  estimates  of 
the  potential-disturbance  I  /e  numbers  and  rigidities  (or,  to  be  more  precise,  their 
real  parts)  for  the  Earth  anc  he  Moon.  The  estimates  are  mutually  consistent.  This 
consistency  unravels  when  ii  comes  to  estimating  the  imaginary  parts  of  k  and  /r. 

The  principal  effect  of  '^[k]  in  the  theory  of  the  libration  of  the  Moon  is  that 
the  node  of  the  lunar  equator  on  the  ecliptic  advances  for  positive  ‘^[k].  Let  SCI 
be  the  offset  in  the  node;  then  I  SCI  =  ^[Ar]  x  208"  [Eckhardt,  1981].  Using  (30), 
Q[fc]  =  -k/<;  =  -0.25  X  10'^,  so,  by  the  nominal  solution,  ISCl  decrea.ses  by 
0".052  leading  to  a  -45  cm  cos  F  term  in  the  latitude  libration.  Also,  Q  =  kq/ k  = 
100  according  to  the  nominal  solution  (1".0  residual  free  libration  and  Eckhardt 
resonance),  but  if  the  Moons  or  Migus  resonance  is  used  instead,  then  Q  a  55. 
Estimates  of  Q  within  the  plausible  range  of  resonances  and  for  different  residuals 
can  be  effected  by  using  Figure  6  which  is  a  plot  of  |7-i(/i)|  as  a  function  of  the 
resonance  period  for  Q  =  30,40, 50, . . . ,  100.  For  low  Q,  as  at  Q  =  30,  rj  ( /i) 
oscillates  about  the  steady  state  solution  (10).  As  an  example,  one  can  see  from  this 
figure  that  with  a  resonance  of  38.658  sidereal  months  and  Q  =  40,  the  residual 
is  jri(/i)l  =  1".15,  and  with  the  same  resonance  and  Q  =  30,  the  residual  is 
|ri(/i)|  =  0".80.  (The  amplitude  in  the  Eckhardt  [1982]  table  is  0".274.)  Thus  a  Q 
as  low  as  40,  and  perhaps  even  30,  may  not  be  inconsistent  with  the  observed  free 
libration. 

Williams,  Newhall  and  Dickey  [1987]  estimate  Q  for  a  monthly  period  to  be 
approximately  30;  in  effect,  they  detected  and  isolated  a  1 .5  m  cos  F  term  in  the 
latitude  libration.  Seismic  estimates  of  the  lunar  rigidity  Q  are  depth  dependent, 
varying  from  about  3000  in  the  region  between  60  and  400  km  depth  to  1000  or 


LUNAR  FREE  LIBRATION  IN  LONGITUDE 


323 


less  for  depths  exceeding  1 100  ktn  [Goins,  Dainty  and  Toksoz,  1981  ].  The  broad 
range  of  lunar  Q  estimates  is  not  unlike  the  range  of  Q  estimates  for  the  Earth. 
Seismic  rigidity  (shear  wave)  Q's  are  depth  dependent,  varying  from  about  100 
at  the  top  of  the  mantle  to  about  2000  at  the  base  of  the  mantle  [Stacey,  pp.  300- 
307,  1977].  Estimates  of  Q  from  polar  motion  data  range  between  30  and  600 
[Munk  and  MacDonald,  pp.  167-174,  1960;  Stacey,  p.  67,  1977).  The  Earth  has 
a  liquid  core  and  a  more  pronounced  stratification  than  the  Moon,  so  it  is  a  mote 
complicated  body.  It  is,  however,  much  easier  to  make  Q-sensitive  measurements 
of  the  Earth  than  of  the  Moon.  That  the  seismic  Q’s  of  the  Moon  are  much  larger 
than  its  rotational  Q's  is  not  an  inconsistency;  rocks  have  different  viscoelastic 
properties  at  ~  1  Hz  than  at  one  month  or  three  years.  The  rotational  Q's  ate 
roughly  bracketed  somewhere  between  30  and  100;  to  some  geophysicists,  that’s 
pinpoint  precision. 


7.  Conclusion 

The  planetary  perturbation  term  ELP-2000  Sequence  Number  16331  [Chapront  and 
Chapront-Touze  ,  1983]  modulates  the  lunar  longitude  with  an  amplitude  of  only 
0".0021  but,  because  it  passed  through  resonance  33  to  66  thousand  years  ago,  it  is 
the  source  of  a  ~  1"  lunar  libration  term  that  has  been  considered  a  free  libration. 
In  the  sense  that  this  libration  term’s  phase  is  a  free  parameter  unrelated  to  other 
parameters  of  lunar  motion,  this  is  true;  but  in  the  sense  that  the  term’s  amplitude 
is  not  a  free  parameter  but,  in  fact,  a  sensitive  indicator  of  lunar  dissipation  at  the 
38.65  sidereal  month  resonant  period  and  its  excitation  source  is  known,  regarding 
the  term  as  a  free  libration  is  misleading  by  accepted  standards.  The  definition  of 
a  free  libration  is  hereby  broadened  to  include  this  term. 

Three  libration  theories  differ  in  their  resonance  periods  for  the  libration  in 
longitude  by  less  than  0.02%  but  the  amplitude  of  the  free  libration  in  longitude  is 
very  sensitive  to  which  theory  is  used;  and,  of  course,  it  is  sensitive  to  the  imaginary 
part  of  the  potential-disturbance  Love  number  at  resonance.  The  Moon  is  a  fairly 
uniform  body,  and  its  rigidity  and  density  ate  well  constrained  and  so,  therefore, 
is  the  real  part  of  the  Lov^'  number.  The  Q  of  the  lunar  rigidity  is  the  dissipation 
parameter  that,  with  a  choice  of  resonance,  can  be  estimated  from  the  free  libration 
amplitude.  Allowing  for  uncertainties  in  the  amplitude  and  resonance,  the  Moon’s 
Q  still  falls  within  a  region  that  is  eniin  'y  reasonable  by  geophysical  standards. 

Since  an  astronomic  excitation  source  for  the  free  libration  in  longitude  has 
been  found,  there  is  less  justification  to  posit  a  source  within  the  Moon,  e.g.,  a 
turbulent  fluid  core  that  interacts  with  the  mantle  [Yoder,  1981].  Thus,  the  case  for 
a  fluid  lunar  core  is  weakened. 
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